clearvars -except Globaloption option 


load MAT/benchmark_US_before1946_demean.mat;
s_oneVAR_break=s;

load MAT/benchmark_US_before1946.mat;

f=figure;
f.Position=[1000 918 560 300];

plot(1793:1946, s_oneVAR_break,'m--','LineWidth',3)
hold on
plot(1793:1946,s,'red','LineWidth',3)
grid on
load '../before1946_twoVAR/MAT/benchmark_US_pre1860.mat'
spre1860=s;
load '../before1946_twoVAR/MAT/benchmark_US_post1860.mat'
s_twovar=[spre1860 s];
plot(1793:1946,s_twovar, 'g--','LineWidth',3)
hold on
load MAT/benchmark_US_before1946.mat;
plot(1793:1946, gdebt(2:end),'blue', 'LineWidth',3);
hold on
ylim([-.5 1.2])
xlim([1793 1946])
legend('One VAR with Break at 1860','One VAR with no Break','Two VARs with Break at 1860','Debt/GDP','Location','northwest')
xlabel('Year')
ylabel('PV(Surplus)/GDP')

f.PaperSize = [6 3.75];

print('../../results/figures_main/US/pre1946/cs_twoVAR_comparison','-dpng','-r1000');
